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~£T) '. Abstract 

1 This paper considers orthogonal frequency division multiple access (OFDMA) systems with simultaneous 

wireless information and power transfer. We study the resource allocation algorithm design for maximization 
of the energy efficiency of data transmission (bits/Joule delivered to the receivers). In particular, we focus 
on power splitting hybrid receivers which are able to split the received signals into two power streams for 
concurrent information decoding and energy harvesting. Two scenarios are investigated considering different 
power splitting abilities of the receivers. In the first scenario, we assume receivers which can split the received 
power into a continuous set of power streams with arbitrary power splitting ratios. In the second scenario, we 
O ■ examine receivers which can split the received power only into a discrete set of power streams with fixed power 

splitting ratios. In both scenarios, we formulate the corresponding algorithm design as a non-convex optimization 
£> ' problem which takes into account the circuit power consumption, the minimum data rate requirements of delay 

constrained services, the minimum required system data rate, and the minimum amount of power that has 
to be delivered to the receivers. Subsequently, by exploiting fractional programming and dual decomposition, 



o 

£f~) ' suboptimal iterative resource allocation algorithms are proposed to solve the non-convex problems. Simulation 

o r 

_ results illustrate that the proposed iterative resource allocation algorithms approach the optimal solution within 

a small number of iterations and unveil the trade-off between energy efficiency, system capacity, and wireless 
power transfer: (1) wireless power transfer enhances the system energy efficiency by harvesting the energy in the 
j_j ■ radio frequency, especially in the interference limited regime; (2) the presence of multiple receivers is beneficial 



for the system capacity, but is not necessarily beneficial for the system energy efficiency. 
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I. Introduction 

Orthogonal frequency division multiple access (OFDMA) has been widely adopted as the air interface 
in high speed wireless multiuser communication networks, due to its immunity to channel delay 
spread and flexibility in resource allocation. In an OFDMA system, a wideband frequency selective 
spectrum is converted into a number of orthogonal narrowband frequency flat subcarrier channels which 
facilitates the multiplexing of users' data and the exploitation of multiuser diversity. On the other 
hand, next generation communication systems are expected to support multiple users and to guarantee 
quality of service (QoS). The increasing demand for high data rate and ubiquitous services has led 
to a high energy consumption in both transmitter(s) and receiver(s). Unfortunately, portable mobile 
devices are usually only equipped with limited energy supplies (batteries) which creates bottlenecks in 
perpetuating the lifetime of networks. Besides, the battery capacity has improved at a very slow pace 
over the past decades and is unable to satisfy the new energy requirements [3]. Consequently, energy- 
efficient mobile communication system design has become a prominent approach for addressing this 
issue in energy limited networks Hl-G] . Specifically, an enormous number of technologies/methods 
such as energy harvesting and resource allocation optimization have been proposed in the literature 
for improving the energy efficiency (bits-per- Joule) of wireless communication systems. Among the 
proposed technologies, energy harvesting from the environment is particularly appealing as it constitutes 
a perpetual energy source. More importantly, it provides self-sustainability to systems and is virtually 
free of cost. 

In practice, numerous renewable energy sources can be exploited for energy harvesting, including 
solar, tide, geothermal, and wind. However, these natural energy sources are usually location, weather, 
or climate dependent and may not always be available in enclosed/indoor environments or suitable for 
mobile devices. On the other hand, wireless power transfer technology, which enables the receivers to 
scavenge energy from propagating electromagnetic waves (EM) in radio frequency (RF), has gained 
recent attention in both industry and academia |[8ll- [fT4ll . Indeed, RF signals carry both information 
and energy simultaneously. Thus, the RF energy radiated by the transmitter(s) can be recycled at the 
receivers for prolonging the lifetime of networks. Yet, the utilization of EM waves as a carrier for 
simultaneous information and power transfer poses many new research challenges for both resource 
allocation algorithm and receiver design. In IfTOll and IfTTTl . the fundamental trade-off between wireless 
information and wireless power transfer was studied for flat fading and frequency selective fading, 
respectively. Specifically, an ideal receiver was assumed in ifTOl and iPTO such that information decoding 
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and energy harvesting can be performed on the same received signal which is not possible in practice yet. 
As a compromise solution, three different types of the receivers, namely power splitting, separated, 
and time -switching receivers, were proposed in lfT2ll . |fT3l . In particular, the power splitting receiver 
splits the received power into two power streams with a certain power splitting ratio for facilitating 
simultaneous energy harvesting and information decoding in the receiver. The authors in W2\ and 
[TT3l investigated the rate-energy regions for different types of receivers in two-user and point-to- 
point single carrier systems, respectively. However, the problem formulations in iTPTl and lfT3l do 
not take into account the heterogeneous data rate requirements of users and the results may not be 
to guarantee the requirements and applicable to multi-carrier systems with arbitrary numbers of users. 
On the other hand, the authors in lfl4ll focused on the resource allocation algorithm design for a 
point-to-point single user system with power splitting receiver in ergodic fading channels. Yet, the 
assumption of channel ergodicity may not be well justified for delay sensitive services in practice 
since the transmitted data symbols of these services experience slow fading. Besides, the high power 
consumption in electronic circuitries and RF transmission have been overlooked in lfT0l - lfT4l but play 
an important role in designing energy efficient communication systems lPT5Tl . IfToll . In other words, the 
energy efficiency of the systems in these works remains unknown. In addition, a portion of RF energy 
can be harvested by the RF energy harvesting receivers for improving the system energy efficiency. 
Yet, the system models adopted in the literature [TT0l - lfT4ll do not consider the energy recycling process 
from an energy efficiency point of view. By incorporating the circuit power consumption and the RF 
energy harvesting ability of receivers in the problem formulation, we presented in Q]| and |2|] two 
energy efficient resource allocation algorithms for multicarrier systems assuming separated receivers 
and power splitting receivers, respectively. Nevertheless, El and O do not fully exploit the degrees of 
freedom in resource allocation since data multiplexing of different users on different subcarriers was 
not considered. Moreover, the algorithm proposed in [0 incurs a high computation complexity at the 
transmitter since the optimal power splitting ratio is found via a full search over a continuous variable. 
Furthermore, the power splitting ratio may take only discrete levels in practice and the results in [0, 
which were designed for continuous power splitting ratios, are no longer applicable. 

In this paper, we address the above issues. To this end, we formulate the resource allocation algorithm 
design for energy efficient communication in OFDMA systems with simultaneous wireless information 
and power transfer as an optimization problem. In particular, we focus on the algorithm design for power 
splitting receivers and consider both continuous and discrete power splitting ratios. Besides, users data 
multiplexing on different subcarriers is incorporated in the problem formulation. The resulting non- 
convex optimization problems are solved by iterative algorithms which combine nonlinear fractional 
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Fig. 1. An OFDMA downlink communication system with K = 3 mobile receivers. The upper half of the figure illustrates a block 
diagram of the transceiver model for wireless information and power transfer. 

programming and dual decomposition. Simulation results illustrate an interesting trade-off between 
energy efficiency, wireless power transfer, and multiuser diversity. 

II. System Model 

In this section, we present the adopted OFDMA signal models and the model for the hybrid infor- 
mation and energy harvesting receiver. 

A. OFDMA Channel Model 

We consider an OFDMA downlink system in which a transmitter services K mobile receivers. In 
particular, each mobile receiver is able to decode information and harvest energy from the received radio 
signals. All transceivers are equipped with a single antenna, cf. Figure [T] The total system bandwidth is 
B Hertz and there are uf subcarriers. We focus on quasi-static block fading channels and assume that 
the downlink channel gains can be accurately obtained by feedback from the receivers. The downlink 
received symbol at receiver k G {1, . . . , K} on subcarrier i e {1, . . . , rip} is given by 



Yi,k — a/ Pi,kh9kHi^X ik + I ik + Z^ k + 



(1) 



where X ijk , P i>k , and H i k are the transmitted data symbol, the transmitted power, and the multipath 
fading coefficient from the transmitter to receiver k on subcarrier i, respectively. l k and g k represent 
the path loss and shadowing attenuation from the transmitter to receiver k, respectively. Z s ik and Z^ k 
are additive white Gaussian noises (AWGN) originating from signal processing and the antenna on 
subcarrier i of receiver k, respectively. They are modeled as Gaussian random variables with zero 
means and variances a 2 zS and a 2 za , respectively, cf. Figured] h,k is the received aggregate co-channel 
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interference on subcarrier i of receiver k with zero mean and variance erf. which is emitted by 
unintended transmitters sharing the same frequency channel. 

B. Hybrid Information and Energy Harvesting Receiver 

In practice, the model of an energy harvesting receiver depends on its specific implementation. For 
example, electromagnetic induction and electromagnetic radiation are able to transfer wireless power 
ifTTTl . Ifl4l . Nevertheless, the associated hardware circuitries in the receivers and the corresponding energy 
harvesting efficiency can be significantly different. Besides, the signal used for decoding the modulated 
information cannot be used for harvesting energy due to hardware limitation lfl4ll . In order to isolate 
the resource allocation algorithm design from the specific hardware implementation details, we do not 
assume a particular type of energy harvesting receiver. In this paper, we focus on receivers which 
consist of an energy harvesting unit and a conventional signal processing core unit for concurrent 
energy harvesting and information decoding, cf. Figure [TJ In particular, we adopt a receiver which 
splits the received signal into two power streams Ifl4ll in the RF front end with power splitting ratios 3 
p\ k and pf k . Subsequently, the two power streams with power splitting ratio pf k and p\ k are used 
for harvesting energy and decoding the modulated information in the signal, respectively. Indeed, by 
imposing power splitting ratios^ of p\ k = 1, pf k = and p\ k = 0,pf k = 1, the hybrid receiver 
reduces to a tradition information receiver or energy harvesting receiver, respectively. Furthermore, we 
assume that the harvested energy is used to replenish a rechargeable battery at the receiver. Besides, 
each receiver has a fixed power consumption of Pc R Watts which is used for maintaining the routine 
operations in the receiver and is independent of the amount of harvested power. We note that in practice 
the receivers may be powered by more than one energy source and the harvested energy may be used 
as a supplement for supporting the energy consumptions^ of the receivers IfTTll . 

III. Resource Allocation - Continuous Set of Power Splitting Ratios 

In this section, we consider the resource allocation algorithm design for maximizing the system energy 
efficiency for the case of a continuous set of power splitting ratios. The derived solution provides not 
only a useful guideline for choosing a suitable number of discrete power splitting ratios in the power 

'indeed, p\ k and pf k represents the fraction in splitting the received power of user k on subcarrier i for information decoding and 
energy harvesting, respectively. Yet, we follow the convention in the literature fl41 and adopt the term "power splitting ratios" in the 
paper. 

2 In this paper, a perfect passive power splitting unit is assumed; i.e., the power splitting does not incur any power consumption and 
does not introduce any power loss or signal processing noise to the system. 

3 In this paper, the unit of Joule-per-second is used for energy consumption. Thus, the terms "power" and "energy" are interchangeable. 
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splitting unit, but also serves as a performance benchmark for the case of discrete power splitting ratios. 
Now, we define the system energy efficiency by first introducing the weighted system capacity and the 
power dissipation of the system. 

A. System Energy Efficiency 

Assuming the availability of perfect channel state information (CSI) at the receiver, the channel 
capacity between the transmitter and receiver k on subcarrier i with subcarrier bandwidth W = B/uf 
is given by3 



PiJk9k\Hi, 



12 



C ijfc = W \og 2 [I + P iM T hk ) and T i>k = f /TTT'f ' 2 » & 

where Pi^i,k is the received signal-to-interference-plus-noise ratio (SINR) on subcarrier i at receiver 
k. The weighted system capacity is defined as the aggregate number of bits delivered to K receivers 
and is given by 

Tip K 

u{v,s, P ) = j2Y, akSi ' kCi ' k [ bits/s ]> (3) 

i=l k=l 

where V = {Pi, k > 0, Vz, k} is the power allocation policy, S = {s iyk = {0, 1}, Vz, k] is the subcarrier 
allocation policy, and p = {pf k ,pf k > 0,Vz,A;} is the power splitting policy with variables p\ k and 
pf k introduced in Section ITl-B I a k > 0, Vfc, is a non-negative constant which accounts for the priorities 
of different receivers and is specified by the application layer. In practice, proportional fairness and 
max-min fairness can be achieved by varying the values of a k over timej [20]. On the other hand, for 
facilitating the resource allocation algorithm design, we incorporate the total power consumption of the 
system in the optimization objective function. In particular, the power consumption of the considered 
system, Utp(V,S, p), consists of five terms and can be expressed as: 

K n F K K 

U TP (V,S,p) = P CT + KPc R + J2J2 £P i,kSi,k-J2^^-J2^^ [ Joule/s ]' ( 4 ) 

k=l i=l k=l k=l 

np K 

where Q Dk = ^(^PijSi,^hgk\Hi !k \ 2 r)kpf,k ( 5 ) 

i=i j=i 



Power harvested from information signal at receiver k 



up 



and Q h = ^(4>+<>&% • (6) 



i=i 



Power harvested from interference signal and antenna noise at receiver k 



4 The received interference signal Ii on each subcarrier is treated as AWGN in order to simplify the algorithm design 1131 . 1181 
3 Optimizing the value 
to 1191 for further details 



3 Optimizing the value of ak for achieving different system objectives is beyond the scope of this paper. Interested readers may refer 
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The first three terms in ©, i.e., Pq t + KPc R + Ylk=i ^i=i £ Pi,kSi,k, represent the power dissipation 
required for supporting reliable communication. P Ct > is the constant signal processing circuit 
power consumption in the transmitter, caused by filters, frequency synthesizer, etc., and is independent 
of the power radiated by the transmitter. Variables KPc R and J2k=i ^2i=i £p> i,kSi,k denote the total 
circuit power consumption in the K receivers and the power dissipation in the power amplifier of 
the transmitter, respectively. To model the power inefficiency of the power amplifier, we introduce a 
multiplicative constant, e > 1, for the power radiated by the transmitter in © which takes into account 
the joint effect of the drain efficiency and the power amplifier output backoff [fTSTl . For example, if 
e = 10, then 10 Watts of power are consumed in the power amplifier for every 1 Watt of power radiated 
in the RF; the wasted power during the power amplification is dissipated as heat. On the other hand, the 
last two terms in ©, i.e., — J2k=i Qo k — J2k=i Qh> represent the harvested energy at the K receivers. 
The minus sign in front of ^2 k=1 Qo k in © indicates that a portion of the power radiated in the RF 
from the transmitter can be harvested by the K receivers. Besides, < r/k < 1 in © and © is a 
constant which denotes the energy harvesting efficiency of mobile receiver k in converting the received 
radio signal to electrical energy for storage. In fact, the term r]klk9k\Hi,k\ 2 pfk in © can be interpreted 
as a frequency selective power transfer efficiency for transferring power from the transmitter to receiver 
k on subcarrier i. Similarly, the minus sign in front of Ylk=\ Qh m © accounts for the ability of 
receivers to harvest energy from interference signals. The weighted energy efficiency of the considered 
system is defined as the total average number of bits successfully conveyed to the K receivers per Joule 
consumed energy and is given by 

U eff (V,S,p) = r ^' S ' P) , [bits/Joule]. (7) 
UTP{r, o, p) 

It is worth mentioning that unlike in other system models used in the literature lfT8l . [12111 . here the 
unintended interference signal may be beneficial to the system performance as far as the system energy 
efficiency is concerned. Although strong interference impairs the channel capacity, strong interference 
can act as a vital energy source which supplies energy to the receivers and facilities energy savings in 
the system. 

Remark 1: Mathematically, it is possible that Utp(V, S, p) takes a negative value. Yet, Utp(V, S, p) > 
always hold in practical communication systems due to the following reasons. First, it can be observed 

that YJiZi Y!k=i eP i,k s i,k > YJiZiY,k=i p i,k s i,k > X^iQ^fc* where the strict inequality is due to the 
second law of thermodynamics from physics. In particular, the communication channel between the 
transmitter and the K receivers is a passive system which does not introduce extra energy to the signals. 
Besides, the desired signal energy received at the receiver is attenuated due to path loss and energy 
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scavenging inefficiency. Second, the interference power is controlled (via regulation) to a reasonable 
level Il22l for reliable communication which results in P Ct + KP Cr 3> Ylk=\ Qh- 

B. Optimization Problem Formulation 

The optimal power allocation policy, V*, subcarrier allocation policy, S*, and power splitting policy 
p*, can be obtained by solving the following optimization problem: 

max U eff {V,S,p) 

V,S,p 

n p K 

s.t. CI: Q Dk +Qi k > P££ h ,Vk, C2:J2Y, P ^' k ^ P —> 

i=i fe=i 

np K np K 

C3: Pq t + £P%,k s i,k < Ppd C4: Sj,fcCi,fc > -R mm , 

i=i fc=i j=i fc=i 

" F 

C5: ^ Si )fe Ci, fe > R mink ,Vk G £>, C6: P i(fc > 0, Vi, fc, 



i=l 



A' 



C7: 5 i)fc = {0,l},Vi J A;, C8: ^s i)fc <l, Vi, 

fc=i 

C9: pf < pf fc < p§, Vz, fc, CIO: p[ < ^ fc5i)fc < p 7 a , Vi, fc, 

Cll:p^, fc + pf fc <l, Vi,*, C12:pf fe = pf fc , VA;,i^j. (8) 

Variable in CI is a constant which specifies the minimum required power transfer to receiver k. 

The value of P max in C2 puts an upper limit on the power radiated by the transmitter. The value of 
-Pmax is a constant which depends on the hardware limitations of the power amplifier. C3 limits the 
maximum power supplied by the power grid for supporting the power consumption of the transmitter 
to Ppg, cf. Figured] Rm ™ in C4 is the minimum required data rate of the system. Although -R rmn is not 
an optimization variable in this paper, we can strike a balance between the system energy efficiency and 
the total system throughput by varying its value. In particular, when i? m i n is increasing, the resource 
allocator may increase the transmit power for satisfying the higher data rate requirement by sacrificing 
the system energy efficiency. C5 is the minimum required data rate -R m i nfc for the delay constrained 
services of receiver k, which is specified by the application layer, and V denotes a set of receivers 
having delay constrained services. C6 is the non-negative orthant constraint of the power allocation 
variables. C7 and C8 indicate that each subcarrier can be allocated to at most one receiver exclusively in 
conveying information; inter-user interference is avoided in the system. Besides, C8 indicates that some 
subcarriers can be excluded from the subcarrier selection process for energy efficiency maximization. 
Boundary variables p§ and pf in C9 denote the constant upper and lower bounds of the power splitting 
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ratio for harvesting energy, respectively. The bounds are used to account for the limited capability of the 
receivers in splitting the received power. Similarly, p\j and p£ in CIO denote the constant upper and lower 
bounds of the power splitting ratio for information decoding where p§ + p£ = 1 and pf + p{j = 1. CI 1 
is the power splitting constraint of the hybrid information and energy harvesting receivers. Its physical 
meaning is that the power splitting unit, see Figured! is a passive device and no extra power gain can be 
achieved by the power splitting process. C12 constrains the power spitting ratio for energy harvesting, 
pf k , such that it is identical for all subcarriers in each receiver. Theoretically, pf k can be different across 
different subcarriers. However, in this case, an analog adaptive passive frequency selective power splitter 
is required which results in a high system complexity. Therefore, we consider a more practical scenario 
where pf k = p^ fc ,V7c,Vj ^ i. We note that we do not impose p\ k = p^ k ,Vk,Vj ^ i, in the problem 
formulation since different subcarriers can be used for multiplexing the data of different receivers. 

C. Solution of the Optimization Problem 

The key challenge in solving ([8]) is the lack of convexity in the problem formulation. In particular, 
the objective function in ([8]) is a ratio of two non-convex functions which generally results in a non- 
convex function. Besides, constraints C1-C12 do not span a convex solution set due to the integer 
constraint for subcarrier allocation in C7 and the coupled optimization variables in CI. In general, 
there is no standard approach for solving non-convex optimization problems. In the extreme case, an 
exhaustive search or branch-and-bound method is needed to obtain the global optimal solution which 
is computationally infeasible even for a small K and n F . In order to make the problem tractable, we 
transform the objective function and approximate the transformed objective function in order to simplify 
the problem. Subsequently, we use the constraint relaxation approach for handling the integer constraint 
C7 to obtain a close-to-optimal resource allocation algorithm. Next, we introduce the objective function 
transformation via a parametric approach from nonlinear fractional programming. 

D. Transformation of the Objective Function 

For notational simplicity, we define T as the set of feasible solutions^ of the optimization problem 
in © spanned by constraints C1-C12. Without loss of generality, we assume that {V, S,p} e J 7 and 
we denote q* as the maximum energy efficiency of the considered system which can be expressed as 

, U{V*,S*,p*) U(V,S,p) 

q — 1 _ — JXl^x - - ( y ) 

U TP (V*,S*,p*) p,s, p Utp(P,S,pY 

6 We assume that the set is non-empty and compact. 
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Now, we introduce the following important Theorem which is borrowed from nonlinear fractional 
programming Il23ll for solving the optimization problem in ©. 

Theorem 1: The resource allocation policy achieves the maximum energy efficiency q* if and only 

if 

max U(P, S, p) - q*U TP (V, S, p) 
V,S,p 

= U(V*,S*,p*)-q*U TP (V*,S*,p*) = 0, (10) 

for U{V, S,p)>0 and U TP (V, S, p) > 0. 

Proof: Please refer to flTJ Appendix A] for a proof of Theorem 1. 

Theorem \T\ provides a necessary and sufficient condition in describing the optimal resource alloca- 
tion policy. In particular, for an optimization problem with an objective function in fractional form, 
there exists an equivalent optimization problem with an objective function in subtractive form, e.g. 
U(V,S, p) — q*UTp(P,S, p) in the considered case, such that both problem formulations lead to the 
same optimal resource allocation policy. Moreover, the optimal resource allocation policy will enforce 
the equality in (flOl) which provides an indicator for verifying the optimality of the solution. As a result, 
we can focus on the equivalent objective function and design a resource allocation policy for satisfying 
Theorem 1 in the rest of the paper. 

E. Iterative Algorithm for Energy Efficiency Maximization 

In this section, an iterative algorithm (known as the Dinkelbach method^ 11231 ) is proposed for solving 
© with an equivalent objective function such that the obtained solution satisfies the conditions stated 
in Theorem 1. The proposed algorithm is summarized in Table |U (on the next page) and the convergence 
to the optimal energy efficiency is guaranteed if the inner problem (fTTT) is solved in each iteration. 

Proof: Please refer to B71 Appendix B] for a proof of convergence. 

As shown in Table HI we solve the following optimization problem for a given parameter q in each 
iteration in the main loop, i.e., lines 3-12: 

max U(V,S, p) — qU T p(V,S, p) 

T,S,p 

s.t. C1-C12. (11) 

We note that for any value of q generated by Algorithm I in each iteration, U (V, S, p) —qUxpiV, S, p) > 
is always valid; negative energy efficiencies will not occur. Please refer to B71 Proposition 3] for a proof. 
In fact, the transformed objective function, i.e., U(V, S, p) — qUrp(V, S, p), has an interesting pricing 

7 We note that the Dinkelbach method is an application of Newton's method for root finding, please refer to 1151 . 1241 for details. 
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TABLE I 



Iterative Resource Allocation Algorithm. 



Algorithm Iterative Resource Allocation Algorithm [Dinkelbach method] 
l: Initialization: L max = the maximum number of iterations and A = the maximum tolerance 

2: Set q = and iteration index j — 

3: repeat {Iteration Process: Main Loop} 

4: For a given q, obtain an intermediate resource allocation policy {V', S', p'} by solving the problem 
in CLD 

5: if U(V, S', pf) - qU TP (V', S', p') < A then 
6: Convergence = true 

7: return {V\S\ p*} = {V, S', p'} and q* = £%;f/} } 
8: else 

9: Set g = ul^ pw) and j=J + 1 
10: Convergence = false 

ii: end if 



12: until Convergence = true or j = L r 



interpretation from the field of economics. In particular, U(V,S, p) indicates the system profit due to 
information transmission while qUTp(V, S, p) represents the associated cost due to energy consumption. 
Besides, the terms Q Dk and Q Ik in U T p(V, S, p) are the corresponding rebate and discount to the energy 
cost via energy harvesting. The optimal value of q represents a scaling factor for balancing the profit 
and cost. 

1 ) Solution of the Main Loop Problem 4771) : The transformed problem has an objective function 
in subtractive form and is parameterized by variable q. Unfortunately, there are still two obstacles in 
tackling the problem. First, the power splitting variables for information decoding and energy harvesting, 
i.e., p\ k and pf k , are coupled with the power allocation variables in both the objective function and 
constraint CI which complicates the solution. Second, the combinatorial constraint C7 on the subcarrier 
allocation variables creates a disjoint feasible solution set which is a hurdle for solving the problem 
via tools from convex optimization. In order to derive a tractable resource allocation algorithm, we 
approximate the transformed objective function in the following. First, we approximate the weighted 
system capacity as 

np K 

U(V,S,p) = Y,Y. akS *' kWlo & ( 1 + ^ l ,k) &U(V,S,p) (12) 

i=l k=l 
np K 

where U(V,S,p) = Yl a ^ W lo §2 (P^k) (13) 

i=l k=l 
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which is a tight approximation for high SINR, i.e., Pi ; k^i,k 3> 1- On the other hand, we adopt a lower 
bound on Utp(V,S, p) in the transformed objective function in (PTTT) : 

np K 

U TP {V, S, p) > U TP (V, S, p) = P Ct + KP Cr + eP ^i,k 

i=l k=l 

K n F K K 

- J2J2{J2 p ^ s ^) lk ^\ H iM 2 vk-J2Q^ (14) 

k=l i=l j=l k=l 

We note that Utp(V,S, p) is obtained by setting pf k = 1 in Utp{V,S, p). Indeed, Utp(V ',<S, p) can 
be interpreted as the use of a theoretical receiver which is able to fully recycle and harvest the energy 
of the signal used for information decoding^. As a result, the transformed objective function can be 
approximated by 

U(V, S, p) - qU TP (V, S, p) < U(V, S, p) - qU TP (V, S, p). (15) 

On the other hand, by exploiting constraint CI 2, we can rewrite constraint CI in the following form 
to remove the associated non-convexity: 

np K np p^cq 

CI: ^ (^/Vr-;.,) //,///, /A.,,'"///, + X>f ifc % > ^PSVA:, (16) 
i=i j=i i=i ' Pl > k 

where the right hand side of the inequality in (PT6l) is due to constraint C12: pf k = pf k = . . . = pf k = 

. . . = pf lp k . Besides, we can further simplify the algorithm design by replacing C12 with the following 

equivalent constraint: 

C13: pf >k = p% V*, r = {2, . . . , n F }. (17) 

Next, we handle the combinatorial constraint in C7 by time-sharing relaxation. In particular, we relax 
the subcarrier selection variable to be a real value between zero and one instead of a Boolean, i.e., 
C7: < s i: k < 1. As a result, s ijk can be interpreted as a time-sharing factor in allocating subcarrier i to 
K receivers for delivering information. In addition, we introduce two new auxiliary variables and define 
them as P^k = Pi,kSi,k and p\ k = p( k Si^- They represent the actual transmitted power and the power 
splitting ratio for information decoding on subcarrier i for receiver k under the time-sharing condition, 
respectively. On the other hand, we replace Cik in C4 and C5 in (fTTT) by Ci k = W\og 2 (— fjfc), 

3 V s i.k 3 / 

fi,fc = T itk \ pi while C i>k > C i;k . We note that although C i>k > C i>k generally holds, C i>k w C i>k is 
r '> k a i,k 

asymptotically tight in the high SINR regime. Since the feasible solution set of the problem with C^ k in 



8 Note that pf k — 1 is only applied to the received energy harvested from the information signal, but not to the portion harvested from 
the interference signals. 
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C4 and C5 is a subset of the original problem |f[ the solution obtained with the Dinkelbach method for 
the approximated objective function can be used as a suboptimal solution to the original optimization 
problem in ([8]). Nevertheless, it will be shown in the simulation section that the proposed resource 
allocation algorithm achieves a close-to-optimal performance in high SINR. Note that by considering 
the transformed optimization problem with the approximated objective function and relaxed constraint 
C7 in each iteration of the Dinkelbach method, the proposed scheme converges to the optimal solution 
of the problem with the approximated objective function introduced in (fl4l) . 

The problem with approximated objective function is jointly concave with respect to (w.r.t.) all 
optimization variables (cf. Appendix) and it can be further verified that the primal problem satisfies 
Slater's constraint qualification. As a result, strong duality holds and solving the dual problem is 
equivalent to solving the primal problem Il25ll . Motivated by this fact, we solve the primal problem by 
solving its dual problem in the following. 

F. Dual Problem Formulation 

In this subsection, the resource allocation policy is derived via solving the dual problem of CCD with 
approximated objective function. For this purpose, we first need the Lagrangian function of the primal 
problem. The Lagrangian of (PTTI) is given by 

C(w,X,j,^,5,v,n,C,V,S, p) 

np K K np 

=J2 J2( ak + 7 + v k) s i,kCi,k - q(Pc T + KP Cr + E £ ^ k 

i=l k=l k=l i=l 

K np K K n F K 



k=l i=l j=l k=l i=l k=l 

np K np K np K 

~ x ( E E p ^ - p — ) -EE ^ (pU + pf,k - 1) - E c* ( E s ^ - \ 

i=l k=l i=l k=l i=l k=l 

K np 

~E E ^> k { p tk - pf,k) - E v k R ™» 

k=l i=2 keV 
K prcq np K np 

-E^(iP - E (E^Ki^v - (18) 

k=l Pl > k i=l j=l i=l 

'in general, the constraint relaxation used in C7 may result in a superset of the feasible solution set. Yet, it will be shown in the next 
section that the optimal subcarrier allocation policy with respect to the approximated objective function takes values of either zero or 
one on each subcarrier. In other words, the subcarrier allocation policy is a Boolean even though it is allowed to take any value between 
zero and one; i.e., the size of the feasible solution set does not change with the constraint relaxation in C7. 
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where C i)k = Wlog 2 (jr^Lk) and f i)fc = T ijk P i w has elements w k > 0, k G {1,...,A'}, 
and is the Lagrange multiplier vector corresponding to the individual minimum required power transfer 
constraint. A > accounts for the maximum transmit power allowance P max in constraint C2. (3 > 
is the Lagrange multiplier connected to C3 accounting for the power usage from the power grid at 
the transmitter. 7 is the Lagrange multiplier associated with the minimum data rate requirement of 
the system in C4. v denotes the Lagrange multiplier vector connected to the minimum individual 
data rate requirement of receiver k in C5 and has elements v k > 0,k = {1,..., K}. We note that 
Vk = 0, V/c ^ V, for the receivers requiring non-delay constrained services. £ is the Lagrange multiplier 
vector accounting for the subcarrier assignment constraint C8 and has elements Q,i = {1, . . . ,n F }. 
5 is the Lagrange multiplier vector connected to the power splitting constraint Cll and has elements 
$i,k, i — {1) • • ■ ,np}, k — {1, . . . , K}. fi is the Lagrange multiplier vector connected to constraint C12 
and has elements i — {1, . . . , rip}, k — {1, . . . , K}. On the other hand, the boundary constraints 
C6, C7, C9, and CIO on the optimization variables are captured by the Karush-Kuhn-Tucker (KKT) 
conditions when deriving the resource allocation solution in the next section. Thus, the dual problem 
for the primal problem (fTTT) is given by 

x Pi 11 ^ J?& x £ ( w ' A 5 7,/5,^,^,M,C,^,5,p). (19) 

G. Dual Decomposition Solution 

In this section, the optima^ resource allocation policy is obtained via Lagrange dual decomposition. 
Specifically, the dual problem in (fT9l) is decomposed into a hierarchy of two levels. Level 1, the 
inner maximization in (fT9~l) . consists of rip subproblems with identical structure that can be solved in 
parallel. Level 2, the outer minimization in (fT9~l) , is the master problem. The dual problem can be solved 
iteratively. Specifically, in each iteration, the transmitter solves the np subproblems by applying the 
KKT conditions for a fixed set of Lagrange multipliers. Then the solutions of the subproblems are used 
for updating the Lagrange multiplier master problem via the gradient method. 

Level 1 (Subproblem Solution): Using standard convex optimization techniques and the KKT con- 
ditions, for a given q, in each iteration of the Dinkelbach method, the power allocation policy and the 
power splitting policy on subcarrier i for receiver k are given by 

W(a k + -f + v k ) 



Pi.k — s i,kP;. k — S{h 



ln(2) 



Vi, k, (20) 



10 In this section, an optimality refers to the optimality for the problem formulation using the time-sharing assumption and the 
approximated objective function. 
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K 



®i,k = qe + I3e + X- ^(q + w k )l k g k \Hi 



k\ 2 Vk, 



k=l 



E* 
Pl,k 



E* 
Pj,k 



Pi.k 



W k P, 



icq 
minjj 



, V/c, and 



q{(J z a + a Ij JVk - S jtk + fJ. jt k 

s i,kPi,k~ S i,k 



P§ 



Pi 



Wk,j = {2,.. .,n F }, 



2 v /ln(2)(5 4 , fc )(cT 2 2 + cTj 



2(<A 



pi 



(21) 
(22) 

(23) 

(24) 

(25) 

c, if x > 



= y/a%yjm{a k + 7 + v k )(a 2 z a + <r?.J + (5 iik ) M 2 )°z° 

and /ii )fc = 0. Here, operators and [x] c , are defined as [a;] = max{0,x} and [x] 
c, [x\. = x, ifd < x < c, \x\° d = d, if d > x, respectively. The power allocation solution in (l20l) 
is known as multilevel water- filling. In particular, the water-level in allocating power on subcarrier 
i for receiver k, i.e., i^pr^ ' * s not on ^ directly proportional to the priority of receiver k via 
variable a k , but also depends on the channel gains of the other K — 1 receivers via the term Efc=i(<? + 
Wk)lkgk\Hi jk \ 2 r] k . Besides, Lagrange multipliers 7, v k , and w k force the transmitter to transmit with 
a sufficiently high power to fulfill the system data rate requirement, R mm , the individual data rate 
requirements of the receivers having delay constrained services, -R m i nfe , and the minimum power transfer 



requirement, P m ^ k , for receiver k, respectively. Moreover, as can be observed from (1201) . the power 



allocation solution P* k is independent of s iyk which facilitates a simple allocation design. 

On the other hand, the power splitting ratio for information decoding, pf* k , is also in the form 
of water- filling and the water-level depends on the priority of the receiver via a k in (1241) . Besides, 
Lagrange multiplier p ijk affects the power splitting ratio solution for energy harvesting in ([221) and (1231) 
such that pf k will eventually equal p^ k , Vj = {2, . . . ,tif}, as enforced by consensus constraint C12. 
Furthermore, it can be observed from ([241) that if crla+crj, 3> cr 2 s , p\ % — > p T L eventually. If p£ = 0, then 
the solution suggests that when the interference power dominates the power of the signal processing 
noise, an infinitesimally small portion of received power should be used at the receivers for information 
decoding in order to achieve the optimal performance; provided that the data rate constraints C4 and 
C5 are satisfied. In other words, most of the received power at the receivers should be used for energy 
harvesting. This is due to the fact that when o\ a +crj. , S> <y 2 s , the SINR on each subcarrier approaches 



■n , 2 2 — r— -2 t — 2 2 and is independent of p!* k . Besides, U T p(V,S, p) is a monotonically 

decreasing function of pf k . Therefore, pf k — > p§ and p\\ — > p ! L become the optimal power splitting 
policy for energy efficiency maximization. The above fact indicates that in the interference limited 
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regime, a hybrid information and energy harvesting receiver must achieve a higher energy efficiency 
than the traditional pure information receiver. 

On the other hand, subcarrier i is assigned to receiver k when the following selection criterion is 
satisfied: 



1 if k 



arg max 



Mi 



where Mj k 



otherwise 

W (a k + 7 + v k ) log 2 ( P* k l k g k \H iM 
1 



(26) 



log 2 



(27) 



ln(2) m(2)(pf^+a 2 J + a 2 s . 
is the marginal benefit provided to the system when subcarrier i is assigned to serve receiver k. In other 
words, receiver k is selected for information transmission on subcarrier i if it can provide the maximum 
marginal benefit to the system. Besides, if receiver k has a high priority or a stringent individual data 
rate requirement, it will have high values of a k or v k and the resource allocator at the transmitter will 
have a higher preference to serve receiver k with subcarrier i. On the other hand, it can be observed 
from (126*1) that although constraint relaxation is used in constraint C7 for facilitating the design of 
the resource allocation algorithm, the subcarrier allocation policy on each subcarrier for the relaxed 
problem remains a Boolean; time sharing does not occur. 

Remark 2: The above observation for the optimal power splitting policy in the interference limited 
regime is valid for both the problem with the original objective function and the problem with the 
approximated objective function. Indeed, if the case of an interference limited regime is considered, 
the use of the approximated objective function is not necessary. Instead, we can first set pfl — » p§ and 
PiX ~* Pl m tne problem formulation, pf£ and p\\ become constants and the associated non-convexity 
vanishes. Then, we optimize V and S by following a similar approach as used in (f8T) - ([271) . 

Level 2 {Master Problem Solution): The Level 2 master problem in ( fT9l can be solved by using the 
gradient method which leads to the following Lagrange multiplier update equations: 

i=l k=l 

n F K 

0(u+l)= (3{u) - Uu) x (PpG-Pc-J2J2 £p ^)} ' (29) 



A(u + 1)= X(u) x P, 



i=l k=l 



7(u 



np K 



i=l k=l 

tJ „E 



$i,k( u + 1) = kk(u) - U(u) X 1 - p i k - pf k 



,Wi,k, 



(30) 
(31) 
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^r,fc(«+l)= Hr,k{ u ) X (pf,k- Pi,h) j Vr G {2, . . . , np}yk 



W k {u + 1): 



w k {u) - £ 6 («) 



X 



" F 



( Pid) l k9k\Hi, k \ 2 Vk + a l, k ^ 

i=l j=l i=l 
II F 



prcq 
minj, 

Pi,k 



yk, 



v k (u + 1) = v k (u) - £ 7 (u) x f s i>k C ijk - R 



i=l 



yk e X>, 



(32) 



(33) 



(34) 



where index u > is the iteration index and £t(m), t G {1, . . . , 7}, are positive step sizes. Then, the 
updated^ Lagrange multipliers in (|28T)-(|34|) can be used for updating the resource allocation policy 
in (l20l)-(l2~71) via solving the n F subproblems in ( fT9l) . As the primal problem with the approximated 
objective function is jointly concave w.r.t. the optimization variables, it is guaranteed that the primal 
optimal solution can be obtained by solving the problems in Level 1 and Level 2 iteratively, provided 
that the chosen step sizes, £t(ra), are sufficiently small G51l . Il26ll . We note that other than the adopted 
gradient method, different iterative algorithms such as the ellipsoid method can also be used for finding 
the optimal Lagrange multipliers due to the convexity of the dual problem ll25l . 



IV. Resource Allocation Design - Discrete Set of Power Splitting Ratios 

In practice, due to the high complexity associated with a high precision power splitting unit, the 
RF energy harvesting receivers may only be capable of splitting the received power into two power 
streams with a finite discrete set of power splitting ratios. In this section, we design a resource allocation 
algorithm for such receivers. In particular, we assume that there are N distinct power splitting ratios 
for energy harvesting and information decoding at each receiver. Thus, the power splitting ratios for 
energy harvesting and information decoding on subcarrier i for mobile receiver k can be represented 
by the following constraints: 

CU:pf tk = {pf,pf,...,pf",...,pf-}, CIS: fh = {ft,fi,---, (35) 

p^ n and pjj 1 , n E {1,2,..., N}, are the possible power splitting modes for energy harvesting and infor- 
mation decoding adopted in receiver k, respectively. The corresponding resource allocation algorithm 
design can be formulated as the following optimization problem: 



max U eff (V,S,p) 

V,S,p 

s.t. C1-C8, Cll, C12, C14, C15. 



(36) 



"it can be observed that updating is not necessary since it does not affect the result of the subcarrier allocation solution in l !26t . 
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Similar to the case of a continuous set of power splitting ratios, the objective function of the above prob- 
lem formulation inherits the non-convexity of the fractional form of the objective function. Therefore, 
by using Theorem 1, we can transform the objective function from the fractional form to a subtractive 
form and solve the problem via the Dinkelbach method. As shown in Table QJ in each iteration of the 
main loop, i.e., lines 3-12, we solve the following optimization problem for a given parameter q: 

max U(V, S, p) - qU TP (V, S, p) 

V,S,p 

s.t. C1-C8, Cll, C12, C14, C15. (37) 

The additional difficulty in solving the above optimization problem compared to the problem for- 
mulation in ([8]) is the disjoint/discrete nature of the optimization variables pf k and p( k , cf. C14 and 
C15 in (1351) . In general, an exhaustive search is required to obtain the global optimal solution and the 
search space grows in the order of N 2K " F , which may not be computationally feasible for systems of 
moderate size. In the following, we transform (|37|) into an optimization problem with tractable solution 
by exploiting subcarrier time-sharing and the concept of power splitting mode selection. In particular, 
different values of pi" can be treated as different operating modes of receiver k with different equivalent 
SINRs. Then, we combine the subcarrier selection with the operating mode selection by augmenting 
the dimensions of the optimization variables. To this end, we define the channel capacity between the 
transmitter and receiver k on subcarrier % with channel bandwidth W by using power splitting mode n 
as 

c h = Wto & (i + fftr&), rg fc = nIn( ^ l f[^[ , , 08) 

and -P™ fc r^ fc is the received SINR on subcarrier % at receiver k using power splitting mode n for 
information decoding. The weighted system capacity is defined as the total average number of bits 
successfully delivered to the K receivers via the N power splitting modes and is given by 

np K N 

U(V N , ^ = E E E 09) 

i=l k=l n=l 

where V N = {P" fc > 0, Vi, k, n} is the power allocation policy and S N = {s™ k = {0, 1}, Vz, k, n} is the 
subcarrier allocation policy for the case of discrete power splitting ratios. We note that the subcarrier 
allocation policy in (|39l ) incorporates the power splitting mode selection for information decoding. On 
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the other hand, the power consumption of the system can be written as 

n F K N K N K N 

u TP (v N , s n , a n ) = p Ct +kp Cr +J2J2J2 eP iA -EE^-EE^ ( 4 °) 

i=l k=l n=l k=l n=l k=l n=l 

n F K N 

where Q\ = E(EE^X)tel%lV'"< (41) 

i=l j=l m=l 

V v ' 

Power harvested from information signal at receiver k with power splitting mode n 

" F 

and Ql = £(<& + <>Nfc»fe • ^ 



i=i 

v- 



Power harvested from interference and antenna noise at receiver k with power splitting mode n 

Here, = {a% = {0, 1}, V7c,n} is the power splitting ratio selection policy for energy harvesting. In 
the above problem formulation, a% is the optimization variable which captures the selection of the power 
splitting mode p^ n for energy harvesting. Similar to the case of the continuous set of power splitting 
ratios, we consider an approximation of the objective function for facilitating a tractable resource 
allocation algorithm design in the following. First, the system capacity between the transmitter and the 
K mobile receivers can be approximated by 

n F K N 

(a) 



i=l k=l n=l 

Cl k = W\og 2 {p^Tl k ) (43) 

(a) 

and m in (1431) is due to the high SINR assumption. On the other hand, a lower bound for the total 



power consumption of the system is given by 

K rip N 
k=l i=l n=l 

n F K K N K N n F 

- E E ( E E ^) imaw -E E E 

i=l k=l j=l m=l k=l n=l i=l 

where U T p(V N , S N , A N ) is obtained by setting pf"a£ = 1. Then, the objective function for the 
optimization problem with discrete sets of power splitting factors is given by 

U{V N , S N ) - qU TP (V N , S N , A N ) < U(V N , S N ) - qU TP (V N , S N , A N ). (45) 
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A. Optimization Problem Formulation 

With a slight abuse of notation, we reformulate the optimization problem with discrete power splitting 
ratios to be solved in each iteration as follows: 



max U{V N ,S N )- qU TP (V N , S N , A 1 

V N ,S N ,A N 



N 



n F K N 



s.t. CI: J2Qn k + Ql> P ^ k > 

n=l 

n F K N 



C2: EEE^>"^ P -- 

i=l k=l n=l 
n F K N 

C4: E E E S h d l* * Rmin, 



i=l k=l n=l 



C5: 



n F K N 



i=l k=l n=l 



S J\C?, >R mink ,VkeV, 



i=l k=l n=l 



C6:P^>0, Vi,k,n, 



C7: sl k = {0,l},Vi,k,n, 
C9: a£ = {0,l}, Wk,n, 

N 

Cll:p^ fc + ^pf"ar <1, Vi,k,n, C12: ^ pj™ s ™ + p^a n k < 1, Vi,*,n. (46) 



if TV 

fc=l n=l 
TV 

CIO: ^a2 = l,Vfc, 



n=l 
N 



m=l 



m=l 



Constraints C1-C6 have the same physical meanings as in the problem formulation in ([8]). Constraints 
C7-C12 are imposed to guarantee that in each receiver only one power splitting mode can be selected for 
information decoding and energy harvesting. Besides, CI 1 and C12 indicate that no extra power gain can 
be achieved in the power splitting process. On the other hand, the non-convexity of the above problem 
formulation is caused by the coupled optimization variables in CI and the combinatorial constraints in 
C7 and C9. In analogy to the techniques used for solving the optimization problem in (flTT) . we solve 
problem (l46l) in the following two steps. In the first step, we relax constraints C7 and C9 such that 
variables and a™ k can assume any value between zero and one, i.e., C7: < < l,\/i,k,n and 
C9: < a™ k < 1, Vz, k,n. Then s" fc and a™ fc can be interpreted as the time sharing factors for receiver 
k in utilizing subcarrier i with power splitting mode n. We also define a new variable P^ k = s^ k Pp k for 
facilitating the design of the resource allocation algorithm. In fact, P™ k represents the actual transmit 
power of the transmitter for receiver k in subcarrier i if power splitting mode n is used under the time 
sharing condition. In the second step, we replace constraint CI in (l46l) by 



n F 



K N 



II F 



n prcq 



C1 ' : E (E E P ^9k\H hk U + + a]J Vk > V fc, 

Pk 



n. 



i=l j=l m=l 



i=l 



(47) 
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Although CI' is equivalent to CI in (|46|) only if a% takes a binary value, i.e., a£ = {0,1}, and 
J2k=i a k = 1' it be shown in the next section that the optimal solution for a£ has a binary form 
under constraint CI', despite the adopted constraint relaxation. Consequently, the optimization problem 
with the approximated objective function and the constraint relaxation is now jointly concave w.r.t. all 
optimization variablej^. Besides, it satisfies Slater's constraint qualification. Therefore, we can apply 
dual decomposition to solve the primal problem via solving its dual problem. The Lagrangian function 
of the primal problem in (l46l) is given by 



£(w, A, 7 , P, S, v, <p, k, C, V N , S N , A N ) 

np K N K N rip 

=E E E(°* + 7 + u *K*3:* - K Pc - + KPcr + E E E ep ^ 

i=l k=l n=l k=l n=l i=l 

K n F N K K N n F K N 



-E E E ( E pz)'^\H„ t \'n k - E E - - E c ( E E •& 

fc=l 7=1 n=l j'=l fc=l n=l 7=1 fc=l n=l 

n F K N n F K N N 

~ x ( E E E % - p — ) - E E E ^ (p T k s tk + E - \ 

i=l fc=l 77=1 i=l fc=l 77=1 

rip K N N 

-E E E ( E rf m <* + p E k n < - 1) - E ^ 



<-mm fc 

i=l fe=l n=l m=l fceX> 



-EI««(% - E (E E ^tai^r* - X>i + O 

fc=l n=l i = \ j = \ m=l i=l 

K N n F K N 

-E^^E^-^-^^-^ + EEE^)' ( 4g ) 

fc=l n=l 7=1 fc=l n=l 

where A, /3, and 7 are the scalar Lagrange multipliers associated to constraints C2, C3, and C4 in (|46|) . 
respectively. w,v,<p, d, and k; are the Lagrange multiplier vectors for constraints CI', C5, C8, CIO, 
Cll, and C12 which have elements w% > 0, n = {1, . . . , N}, k = {1, . . . , K}, v k > 0, & > 0, % = 
{1, . . . , rip}, (fi, 5™ k , and K™ fc respectively. We note that there is no restriction on the value of ifi since 
it is associated with the equality constraint C10. Thus, the dual problem is given by 

inin max £(w, A, 7, 0, S, v, <p, k, C, V N , S N , A N ). (49) 

w,A,7,/3,<5,^,C,k>0,v5 V N ,S N ,A N 

B. Dual Decomposition Solution 

By using dual decomposition and following a similar approach as in Section IIII-GL the resource 
allocation policy can be obtained via an iterative procedure. For a given set of Lagrange multipliers 

12 The concavity of the above optimization problem can be proved by following a similar approach as in the Appendix for the case of 
continuous power splitting ratios. 
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{w, A, 7, 0, 8, v, ip, k, £}, the power allocation policy, power splitting policy, and subcarrier allocation 
policy for receiver k using power splitting mode n on subcarrier i are given by 

W(a k + 7 + v k ) 



pn* 
r i,k 



n pn* _ 

i,k ~ ° l , k 



ln(2) 



, Vz', k, 



1 if n = argmax Tj! ,Wk, 

b 

otherwise 

Up 



np np N np 

i=l i=l m=l i=l 



Pk 



V>k, 



1 if n, k — argmax M? 6 , Vz, 

c,b 

otherwise 



(50) 

(51) 
(52) 
(53) 



W(a k + 7 + Vk ) [log 2 (P^lk9k\Hi,k\ 2 ) + log 2 (— 



Pk 



^*p£ n -(E<*)*~ Ci (54) 



m=l 



ln(2) ln(2)(p^+^J + a 2 2 s 
and $j fe is defined in (1251) . The power allocation solution in (1501) has a similar multi-level water filling 
interpretation as in ([201) . The difference between ([201) and (1501) is that the power allocation in (1501) is 
performed w.r.t. each power splitting mode. On the other hand, it can be observed from (|51"I) and (1531) 
that the optimal values of a£ and s" fc are binary numbers, although time sharing relaxation is used for 
facilitating the algorithm design. 

Now, since the dual function is differentiable, we can update the set of Lagrange multipliers for a 
given set of V N , S N , A N by using the gradient method. The gradient update equations are given by 



n F K N 

A(«+l)= A(«) - x (p max - E E E P £ 



i=l k=l n=l 

np K N 



f){u + 1) = [£(«) - x(Pp G -P c -^^£ 

i=l k=l n=l 
n F K N 

7 (« + 1) = r 7 («) - e 3 («) x (i? min -EEE s ^ 



i=l fe=l n=l 
N 
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Similar to the case of a continuous power splitting ratios, updating Q and is not necessary since they 
will not affect the power splitting mode selection and subcarrier allocation in (l5TT) and (1531) . respectively. 

V. Results 

In this section, simulation results are presented to demonstrate the energy efficiency and system 
capacity of the proposed resource allocation algorithms. We consider an indoor communication system 
with a maximum service distance of 10 meters. There are K receivers uniformly distributed between 
the reference distance d = 1 meter and the maximum service distance. Each transceiver is assumed 
to have an effective antenna gain of 12 dB. The TGn path loss model 11271 is adopted with a carrier 
center frequency of 470 MHzo. We note that the wavelength of the carrier signal is 0.6 meter which is 
smaller than minimum distance between the transmitter and receivers. Thus, the far-field assumption of 
the channel model in 11271 holds. The operating bandwidth is B = 20 MHz with np = 128 subcarriers 
which results in a subcarrier bandwidth of W = 156 kHz. All receivers are assumed to have the same 
priority a k = 1, Vk, in resource allocation. We assume that the signal processing noise, cr5L, in every 
receiver is due to thermal noise and quantization noise. Specifically, a 12-bit uniform quantizer is used 
for quantizing the information carried on each subcarrier in the fast Fourier transform computation. As 
a result, the quantization noise power and the thermal noise power on each subcarrier are —47 dBm 
and —112 dBm, respectively. In addition, the antenna noise is set to a 2 a = —44 dBm. The multipath 

i,k 

fading coefficients of the transmitter-receiver communication channels are generated as independent 
and identically distributed (i.i.d.) Rician random variables with Rician factor 6 dB. The shadowing of 
all communication links is set to g k — 1, V7c, to account for the line-of-sight communication setting. 
We assume the maximum power supply of the transmitter is P PG = 50 dBm and the minimum data 
rate requirement of the system is R mm = 50 Mbit/s. Unless further specified, we assume that there is 



one receiver requiring delay constrained service with a minimum data rate requirement of R 



10 



Mbit/s. Besides, the minimum required power transfer and the energy harvesting efficiency are set to 



P 



rcq 
min fe 



dBm, Wk, and rj k = 0.8, respectively. Also, we set the static circuit power consumptions 



l3 The 470 MHz frequency band will be used by IEEE 802.11 for the next generation of Wi-Fi systems 1281 



1:45am, Tuesday 19 th March, 2013 



DRAFT 



x 10 



23 



_ 18 

0) 

I 1.6 

CO 

* 1.4 

o 

c 

33 

jg 1.2 

0) 

2> 1 

03 
C 
0) 

1 08 

>, 

o 0.6 



| 0.4 



0.2 






i_ _ Upper bound performance, o 


= -50 dBm 


Proposed algorithm [, <jp = - 


50 dBm 


^ Proposed algorithm II, = 


-50 dBm 


_ A - Upper bound performance, a 


= -40 dBm 


Proposed algorithm I, op = - 


40 dBm 


Proposed algorithm II, = 


-40 dBm 



4 5 6 

Number of iterations 



10 



Fig. 2. Average system energy efficiency (bits/joule) versus number of iterations with different levels of interference power, 07 , and 
different values of maximum transmit power allowance, P max . There are K = 3 receivers in the system. The dashed solid lines represent 
the upper bound of energy efficiency after algorithm convergence. 

of the transmitter and each receiveo to Pc T = 30 dBm and Pc R = 10 dBm, respectively. Furthermore, 
we assume a class A/B power amplifier with a power efficiency of 16% is used at the transmitter li29l . 
i.e., e = jj^g = 6.25. The average energy efficiency of the system is computed according to © and 
averaged over multipath fading and path loss. In the sequel, the total number of iterations is defined 
as the number of main loops in the Dinkelbach method. For the case of continuous power splitting 



ratios, we set pf} 



1 and pf 



Pi 



0. Besides, there are five power splitting ratios for 



the resource allocation with the discrete set of power splitting ratios: pf k = {1,0.75,0.5,0.25,0} and 
p\ k = {0, 0.25, 0.5, 0.75, 1}. Moreover, the step sizes adopted in (T28T)-(l34l) and (l55l)-(l6TT) are optimized 
for obtaining fast convergence. Note that if the transmitter is unable to meet the minimum required 
system data rate R m \ n , the minimum required individual data rate -R m i nfc , or the minimum required 
power transfer P™^ , we set the energy efficiency and the system capacity for that channel realization 
to zero to account for the corresponding failure. For the sake of illustration, the performance curves 
of the proposed algorithms for the continuous and discrete sets of power splitting ratios are labeled as 
"Proposed algorithm I" and "Proposed algorithm II" in Figures [2]-[6l 

A. Convergence and Optimality of Iterative Algorithm 

Figure [2] depicts the average system energy efficiency of the proposed iterative algorithms for 
different levels of received interference versus the number of iterations. Specifically, we are inter- 
ested in the energy efficiency and convergence speed of the proposed algorithms. We plot the upper 

14 The values of the circuit power consumption used in this paper are for illustration purpose. In practice, the value of circuit power 
consumption depends on the specific hardware implementation and service application. 
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bound performance^ of the system after convergence to illustrate the sub-optimality of the proposed 
algorithms. The results in Figure [2] were obtained by averaging the energy efficiency of the system 
over 100000 independent realizations of multipath fading and path loss attenuation. The dashed lines 
refer to the average energy efficiency upper bound for each case study. After only 5 iterations, the 
iterative algorithms achieve over 95% of the upper bound value for all considered scenarios. Besides, 
the convergence speed of the proposed algorithms is invariant to the interference levels, cr| ifc , which is 
desirable for practical implementation. 

In the following case studies, the number of iterations is set to 5 for illustrating the performance of 
the proposed algorithms. 

B. Energy Efficiency versus Maximum Allowed Transmit Power 

Figure [3] shows the average system energy efficiency versus the maximum transmit power allowance, 
-Pmax, for different received levels of interference, erf. k . It can be observed that the average system energy 
efficiency of the proposed algorithms is a monotonically non-decreasing function of -P max . hi particular, 
starting from a small value of P max , the energy efficiency first quickly increases with an increasing 
P max and then saturates when P max > 18 dBm. This is due to the fact that the two proposed algorithms 
strike a balance between the system ene rgy efficiency and the power consumption of the system. In 
fact, once the maximum energy efficiencyLj of the system is achieved, a further increase in the transmit 
power would result in a degradation in energy efficiency. As expected, proposed algorithm I outperforms 
proposed algorithm II in all cases since the latter algorithm is designed based on discrete sets of power 
splitting ratios which span a smaller feasible solution set compared to algorithm I. Besides, Figure [3] 
reveals that although interference signals can act as a viable energy source to the system, cf. © and 
(Q, strong interference impairs the energy efficiency of the system; the energy harvesting gain due to 
strong interference^ is unable to compensate the corresponding capacity loss. For comparison, Figure 
[3] also contains the energy efficiency of two baseline resource allocation schemes which maximize the 
system capacity (bits/s) under constraints C1-C12 in © for different settings of power splitting ratios. 
Specifically, the baseline I and baseline II algorithms maximize the system capacity w.r.t. {V, S}, but 
for a random value of pf k and a fixed value pf k = 0.5, respectively. Figure [3] shows that the proposed 

15 The upper bound system performance is obtained by directly evaluating the upper bound objective function via replacing the energy 
consumption function Utp(V, S, p) in dl4b by U(V, S, p) for the continuous set of power splitting ratios. We note that the approximation 
in J 12l i is asymptotically tight for high SINR. 

16 The maximum energy efficiency refers to the "maximum" w.r.t. the corresponding problem formulation. 

"Nevertheless, we would like to emphasize that the use of hybrid information and energy harvesting receivers provides a better energy 
efficiency to the system compared to pure information receivers, as discussed before equation i26l . 
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Fig. 3. Average system energy efficiency (bits-per- Joule) versus maximum transmit power allowance, P ma x, for K — 3 receivers and 
different levels of interference power, erf . . 

algorithms I and II outperform baseline algorithms I and II for all considered scenarios. The superior 
performance of the proposed algorithms is attributed to the optimization of the energy efficiency and 
the power splitting ratios. It can be observed that the energy efficiency of the two baseline schemes 
reduces dramatically in the high transmit power allowance regime. The reason is that the baseline 
schemes allow the transmitter to use an exceedingly large power for capacity maximization which has 
a negative impact on the energy efficiency. On the other hand, the performance gains achieved by the 
two proposed algorithms over the two baseline schemes is small when the interference power is high 
and the transmit power is low. This is because the system capacity improvement due to an optimized 
value of p\ k is saturated in the high interference regime, although part of the energy in the RF is 
harvested by the receivers. 

C. Average System Capacity versus Maximum Allowed Transmit Power 

In Figure HJ we plot the average system capacity versus the maximum transmit power allowance, 
-P m ax> for different levels of interference power, aj, k . We compare the two proposed algorithms again 
with the two aforementioned baseline schemes. For P max < 18 dBm, it can be observed that the average 
system capacities of the two proposed algorithms scale with the maximum transmit power allowance 
-Pmax- Yet, the system capacity gain due to a larger P max begins to saturate in the high transmit power 
allowance regime, i.e., P max > 18 dBm. Indeed, the proposed algorithms do not further increase the 
transmit power in the RF when the system capacity gain due to a higher transmit power cannot neutralize 
the associated energy consumption required for boosting the transmit power. On the other hand, the 
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Fig. 4. Average system capacity (bits-per-second) versus maximum transmit power allowance, P ma x, for different power levels of 
interference, a'j. , and K = 3 receivers. 

proposed schemes achieve a superior average system capacity compared to the baseline schemes in the 
low transmit power allowance regimes. These results suggest that the optimization of the power splitting 
ratio plays also a key role in maximizing the system capacity. Nevertheless, the system capacities of 
the two baseline schemes are larger than those of the proposed algorithm in the high transmit power 
regime. This is attributed to the fact that in order to maximize the system capacity for the baseline 
algorithms, the transmitters radiate all the power at every time instant whenever it is available. Yet, 
the better system capacity comes at the expense of a low system energy efficiency, see Figure [3] In 
addition, all the considered algorithms are able to fulfill the system data rate and individual data rate 
requirements in C4 and C5 on average. 

D. Average Harvested Power versus Maximum Allowed Transmit Power 

Figure [5] depicts the average harvested power of the proposed algorithm II versus the maximum 
allowed transmit power, P max , for different levels of interference power, aj . It can be seen that in 
the high transmit power regime, the amount of average harvested power in all considered scenarios is 
saturated. This is because for the proposed algorithm II, the transmitter stops to increase the transmit 
power for energy efficiency maximization. Meanwhile, the average interference power level remains 
unchanged and no extra energy can be harvested in the K receivers. On the other hand, a higher amount 
of power is harvested by the receivers in the system when the interference power levels increases. The 
reason behind this is twofold. First, the increases in interference power level provide some extra energy 
to the system for potential energy harvesting. Second, the strong interference tends to saturate the 
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Fig. 5. Average harvested power (dBm) of proposed algorithm II versus the maximum transmit power allowance, P max , for different 
interference power levels and K = 3 receivers. The double-sided arrow indicates the power harvesting gain due to an increasing 
interference power level. 



SINR on each subcarrier such that it is independent of p(* k , i.e. 



Pi,kPl* k lk9k\Hi,k\ 2 



y Pi,khgk\Hi,k\ 2 



a 2 z a + erf. 3> a 2 s . Thus, using more of the received power for information decoding does not provide a 
significant gain in channel capacity. Consequently, more received power is used for energy harvesting 
to reduce the total energy consumption of the system which enhances the system energy efficiency. 
E. Average Energy Efficiency and System Capacity versus Number of Receivers 

Figure |6(a)| and Figure |6(b)| illustrate the average system capacity and the average system energy 
efficiency of the proposed algorithm II versus the number of mobiles receivers for different interference 
power levels, crf ifc , and different receiver circuit power consumptions, Pc R . In Figure |6(b~)| it can be 
observed that the average system capacity increases with the number of receivers in the system since 
the proposed algorithm is able to exploit multiuser diversity. Specifically, the transmitter has a higher 
chance of selecting a receiver with good channel conditions when more receivers are in the system, 
which results in a system capacity gain. In addition, although a higher interference power level impairs 
the average system capacity, it does not decrease the performance gain due to multiuser diversity as 
can be concluded from the slopes of the curves. Besides, a higher circuit power consumption in the 



receivers does not have a large impact on the average system capacity. On the contrary, Figure |6(a) 
shows that the average system energy efficiency does not necessarily monotonically increase/decrease 
w.r.t. the number of receivers. In fact, for a moderate value of receiver circuit power consumption, e.g. 
Pc R > 5 dBm, for the considered system setting, the energy efficiency of the system first increases and 
then decreases with an increasing number of receivers. The enhancement of energy efficiency is mainly 
due to the multiuser diversity gain in the channel capacity when having multiple receivers. Besides, if 
more receivers participate in the energy harvesting process, a larger portion of energy can be harvested 
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Fig. 6. Average energy efficiency (bits-per- Joule) and average system capacity (bit-per-second) versus the number of receivers K for 
fmax = 18 dBm, different receiver circuit power consumptions, Pc R , and different interference power levels a'j. . 

from the RF signals. Nevertheless, an extra circuit energy consumption is incurred by each additional 
receiver. Indeed, when the number of receivers in the system or P Cr are large, the system performance 
gain due to multiuser diversity is unable to compensate the total energy cost of the receivers since 
KP Cr increases linearly w.r.t. the number of receivers. Hence, the energy efficiency of the system 
decreases with the number of receivers. As a matter of fact, the energy efficiency gain due to the 
additional receivers depends on the trade-off between multiuser diversity gain, the amount of harvested 
energy, and the associated cost in having multiple receivers in the system. In the extreme case, the 
energy efficiency will monotonically increase w.r.t. the number of receivers if P Cr —> 0, provided that 
the optimization problems in © and (l46l) are feasible. 

VI. Conclusions 

In this paper, the resource allocation algorithm design for simultaneous wireless information and 
power transfer in OFDMA systems was studied. We focused on power splitting receivers which are 
able to split the received signals into two power streams for concurrent information decoding and 
energy harvesting. The algorithm design was formulated as a non-convex optimization problem which 
took into account a minimum system data rate requirement, minimum individual data rate requirements 
of the receivers, a minimum required power transfer, and the total system power dissipation. We first 
focused on receivers with continuous sets of power splitting ratios and proposed a resource allocation 
algorithm. The derived solution served as a building block for the design of a suboptimal resource 
allocation algorithm for receivers with discrete sets of power splitting ratios. Simulation results showed 
the excellent performance of the two proposed suboptimal algorithms and also unveiled the trade-off 
between energy efficiency, system capacity, and wireless power transfer. 
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APPENDIX 

A. Proof of Concavity of the Transformed Problem with Objective Function Approximation 

The concavity of the optimization problem with approximated objective function can be proved by 
the following few steps. First, we consider the concavity of function U(V,S,p) on a per subcarrier 
basis w.r.t. the optimization variables P itk , p( k , and pf k . For notational simplicity, we define a vector 

Xi, fc = [Pi, k Plk pfk\ and a function f ijk (^i,k) = Wa k \og 2 f /^'+^^+l 2 ) which takes vector x ^,fc 

\ "i, k ' z a 1% k 2 s / 

as input. Then, we denote by H(/j fc(xj fc)), and n, Ti, and r 3 the Hessian matrix of function fi,k(*-i,k) 
and the eigenvalues of H(/ i)fc (xj ifc )), respectively. The Hessian matrix of function /^(x^) is given by 



H(/i,fc(x i)fc )) 







Q V ' 'i,kJ Q 





(62) 



and the corresponding eigenvalues are T\ = (f> 7 2 , r 2 = — l ' kJ 2 , and r 3 = 0. 

( " fe) ( ) (pf,J 2 M2)(^+pf, fc ^ 



T 2 

Since r a < 0, a = {1,2,3}, therefore H(/ i j c (x i fc)) is a negative semi-definite matrix. In other words, 
function /^(x^) is jointly concave w.r.t. P ijfc , pf fc , and pf k . Then, we can perform the perspective 
transformation on /^(x^) which is given by M ijfc (xj fc ) = s it kfi,k( x i,k/ s i,k)- We note that the perspective 
transformation preserves the concavity of the function ll25ll and Wj^Xj^) is jointly concave w.r.t. P i>k , 
Pi,k> Pi,k> and Subsequently, [/ ("P, S, p) = Y%=i J2k=i Wa k u i,k(*i,k) can be constructed as a non- 
negative weighted sum of u itk (~x i:k ) which guarantees the concavity of the resulting function. Besides, 
Utp(V,S, p) is an affine function of the optimization variables. Thus, U(V,S, p) — qUrpiV \S, p) is 
jointly concave w.r.t. the optimization variables. On the other hand, constraints Cl-Cll (with relaxed 
constraint C7), C14, and C15 span a convex feasible set. As a result, the transformed problem with the 
approximated objective function is a concave maximization problem. 

References 

[1] D. W. K. Ng, E. Lo, and R. Schober, "Energy-Efficient Resource Allocation in Multiuser OFDM Systems with Wireless 
Information and Power Transfer," in Proc. IEEE Wireless Commun. and Networking Conf., 2013. [Online]. Available: 
|http://arxiv.org/abs/1212.3638| 

[2] , "Energy-Efficient Power Allocation in OFDM Systems with Wireless Information and Power Transfer," in Proc. IEEE 

Intern. Commun. Conf., 2013. [Online]. Available: http://arxiv.org/abs/1301.7657 
[3] K. Pentikousis, "In Search of Energy-Efficient Mobile Networking," IEEE Commun. Magazine, vol. 48, pp. 95 -103, Jan. 2010. 
[4] T. Chen, Y. Yang, H. Zhang, H. Kim, and K. Horneman, "Network Energy Saving Technologies for Green Wireless Access 

Networks," IEEE Wireless Commun., vol. 18, pp. 30-38, Oct. 2011. 
[5] O. Arnold, F. Richter, G. Fettweis, and O. Blume, "Power Consumption Modeling of Different Base Station Types in Heterogeneous 

Cellular Networks," in Proc. Future Network and Mobile Summit, 2010, pp. 1-8. 



1:45am, Tuesday 19 th March, 2013 DRAFT 



30 



[6] J. Yang and S. Ulukus, "Optimal Packet Scheduling in an Energy Harvesting Communication System," IEEE Trans. Commun., 
vol. 60, pp. 220-230, Jan. 2012. 

[7] D. W. K. Ng, E. Lo, and R. Schober, "Energy-Efficient Resource Allocation in OFDMA Systems with Large Numbers of Base 

Station Antennas," IEEE Trans. Wireless Commun., vol. 11, pp. 3292 -3304, Sep. 2012. 
[8] F. Zhang, S. Hackworth, X. Liu, H. Chen, R. Sclabassi, and M. Sun, "Wireless Energy Transfer Platform for Medical Sensors and 

Implantable Devices," in Annual Intern, Conf. of the IEEE Eng. in Med. and Biol. Soc, Sep. 2009, pp. 1045-1048. 
[9] V. Chawla and D. S. Ha, "An Overview of Passive RFID," IEEE Commun. Magazine, vol. 45, pp. 11-17, Sep. 2007. 
[10] L. Varshney, "Transporting Information and Energy Simultaneously," in Proc. IEEE Intern. Sympos. on Inf. Theory, Jul. 2008, pp. 

1612 -1616. 

[11] P. Grover and A. Sahai, "Shannon Meets Tesla: Wireless Information and Power Transfer," in Proc. IEEE Intern. Sympos. on Inf. 

Theory, Jun. 2010, pp. 2363 -2367. 
[12] R. Zhang and C. K. Ho, "MIMO Broadcasting for Simultaneous Wireless Information and Power Transfer," in Proc. IEEE Global 

Telecommun. Conf, Dec. 2011, pp. 1 -5. 
[13] L. Liu, R. Zhang, and K.-C. Chua, "Wireless Information Transfer with Opportunistic Energy Harvesting," IEEE Trans. Wireless 

Commun., vol. 3, pp. 345-362, Jan. 2013. 
[14] X. Zhou, R. Zhang, and C. K. Ho, "Wireless Information and Power Transfer: Architecture Design and Rate-Energy Tradeoff," 

submitted for possible journal publication, 2012. [Online], Available: http://arxiv.org/abs/1205.0618 
[15] C. Isheden, Z. Chong, E. Jorswieck, and G. Fettweis, "Framework for Link-Level Energy Efficiency Optimization with Informed 

Transmitter," IEEE Trans. Wireless Commun., vol. 11, pp. 2946 -2957, Aug. 2012. 
[16] G. Miao, N. Himayat, and G. Li, "Energy-Efficient Link Adaptation in Frequency-Selective Channels," IEEE Trans. Commun., 

vol. 58, pp. 545-554, Feb. 2010. 

[17] F. Philipp, P. Zhao, F. Samman, M. Glesner, K. Dassanayake, S. Maheswararajah, and S. Halgamuge, "Adaptive Wireless Sensor 

Networks Powered by Hybrid Energy Harvesting for Environmental Monitoring," in IEEE Intern. Conf. on Inf. and Autom. for 

Sustainability, Sep. 2012, pp. 285 -289. 
[18] A. Iyer, C. Rosenberg, and A. Karnik, "What is the Right Model for Wireless Channel Interference?" IEEE Trans. Wireless Commun., 

vol. 8, pp. 2662-2671, May 2009. 
[19] Y. Ma, A. Leith, M.-S. Alouini, and X. Shen, "Weighted- SNR-Based Fair Scheduling for Uplink OFDMA," in Proc. IEEE Global 

Telecommun. Conf, Dec. 2009, pp. 1-6. 
[20] G. Song and Y. Li, "Cross-layer Optimization for OFDM Wireless Networks-Part II: Algorithm Development," IEEE Trans. Wireless 

Commun., vol. 4, pp. 625-634, Mar. 2005. 
[21] D. Ng and R. Schober, "Resource Allocation and Scheduling in Multi-Cell OFDMA Systems with Decode-and-Forward Relaying," 

IEEE Trans. Wireless Commun., vol. 10, pp. 2246-2258, Jul. 2011. 
[22] "Report of the Spectrum Ffficiency Working," FCC Spectrum Policy Task Force, Tech. Rep., Nov 2002, [Online] 

http://www.fcc.gov/sptf/reports.html 
[23] W. Dinkelbach, "On Nonlinear Fractional Programming," Management Science, vol. 13, pp. 492-498, Mar. 1967. [Online]. 

Available: http://www.jstor.org/stable/2627691 
[24] T. Ibarako, "Parametric Approaches to Fractional Programs," Mathematical Programming, vol. 26, pp. 288 -300, Jan. 1983. 
[25] S. Boyd and L. Vandenberghe, Convex Optimization. Cambridge University Press, 2004. 

[26] S. Boyd, L. Xiao, and A. Mutapcic, "Subgradient Methods," Notes for EE392o Stanford University Autumn, 2003-2004. 
[27] IEEE P802.ll Wireless LANs, "TGn Channel Models", IEEE 802.1 l-03/940r4, Tech. Rep., May 2004. 

[28] H.-S. Chen and W. Gao, "MAC and PHY Proposal for 802.11af," Tech. Rep., Feb., [Online] 
|https://mentor.ieee.org/802. 1 l/dcn/10/1 1- 10-0258-00-00af-mac-and-phy-proposal-for-802- 1 lafipdfl 

[29] F. Wang, A. Ojo, D. Kimball, P. Asbeck, and L. Larson, "Envelope Tracking Power Amplifier with Pre-Distortion Linearization for 
WLAN 802.1 lg," in IEEE MTT-S Intern. Microw. Symposium Digest, Jun. 2004, pp. 1543 - 1546. 



1:45am, Tuesday 19 th March, 2013 



DRAFT 



